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Abstract 

We consider the asymptotic behavior of perturbations of Lax and 
overcompressive type viscous shock profiles arising in systems of regu¬ 
larized conservation laws with strictly parabolic viscosity, and also in 
systems of conservation laws with partially parabolic regularizations 
such as arise in the case of the compressible Navier-Stokes equations 
and in the equations of magnetohydrodynamics. Under the necessary 
conditions of spectral and hyperbolic stability, together with transver- 
sality of the connecting profile, we establish detailed pointwise esti¬ 
mates on perturbations from a sum of the viscous shock profile under 
consideration and a family of diffusion waves which propagate pertur¬ 
bation signals along outgoing characteristics. Our approach combines 
the recent L^’-space analysis of Raoofi |Baj with a straightforward 
bootstrapping argument that relies on a refined description of non¬ 
linear signal interactions, which we develop through convolution es¬ 
timates involving Green’s functions for the linear evolutionary PDE 
that arises upon linearization of the regularized conservation law about 
the distinguished profile. Our estimates are similar to, though slightly 
weaker than, those developed by Liu in his landmark result on the case 
of weak Lax type profiles arising in the case of identity viscosity !Liu3j. 
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1 Introduction 


Consider a “viscous shock profile”, or traveling-wave solution 
(1.1) u{x,t) = u{x — st), lim u{z) = u±, 

2—>±oo 


of a second-order system of conservation laws 


( 1 . 2 ) 


ut + F{u)^ = 

X eR-,u, F elSF; B e 


for which we take 


(1.3) 







G M"" ”, G M”, and 


(1.4) Reab2 > 0, 

with 9 > 0 and where a denotes spectrum. 

Equations of form (HI arise in a wide variety of applications, includ¬ 
ing such well known examples as the compressible Navier-Stokes equations, 
the equations of magnetohydrodynamics (see |Z.4j h and the equations of 
viscoelastic materials. The long time behavior of solutions of (ll.2|l is of¬ 
ten characterized by a family of viscous shock prohles (EH), which serves as 
a regularization for the shock pattern solution of the associated hyperbolic 
problem 


( 1 . 6 ) 


Ut + F{u)tt = 0. 


In this context, we expect that a viscous prohle will persist in the pattern 
only if it is individually stable to small fluctuations in initial conditions, and 
consequently the stability of viscous shock prohles has long been a subject of 
considerable interest and research effort (see [ZH] and the references therein). 

An important advance in the study of such waves and their stability was 
the identihcation in [ZHj of a spectral criterion based on the Evans function 
(see particularly [Arm H EH m EHl IMaZ.ai ETX ESI [ZH]). Srlehy, the 
Evans function, typically denoted D{X), serves as a characteristic function 
for the linear operator L that arises upon linearization of (EH about the 
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stationary profile u{x). More precisely, away from essential spectrum, zeros 
of the Evans function correspond in location and multiplicity with eigenvalues 
of L |A(hTl IGZl IZHj . It was shown in [ZH] and jMaZ.dj . respectively for the 
strictly parabolic and real viscosity cases, that linearized orbital 

stability of the prohle, p > 1, is equivalent to the Evans function condition, 

(P) There exist precisely i zeroes of D{-) in the nonstable half-plane 
3?A > 0, necessarily at the origin A = 0. 

Here, £ is the dimension of the manifold connecting U- and u+. 

Stability criterion {T>) has been shown to hold in all cases for small am¬ 
plitude Lax shocks arising in strictly parabolic systems [EiEl 10331 m 
IKMl IKMNl IMNl [FT] . as well as for large-amplitude shocks in such cases 
as Lax type waves arising in isentropic Navier-Stokes equations for the 
gamma-law gas as 7 —1 |MNj . and undercompressive shocks arising in 
Slemrod’s model for van der Waal gas dynamics |Z.4j (see [S] for Slemrod’s 


model). More generally, condition (T>) can be verified by numerical calcu¬ 
lation inn ID inzi innoi. In the case of strictly parabolic systems such 
spectral stability, along with standard technical hypotheses on F and B (see 
(7i0)-(7i4), in Remark ll. 71 below. 1 has been shown sufficient for establishing 
nonlinear stability of Lax, under-, over-, and mixed under-overcompressive 
shock prohles [HZj . while for mixed hyperbolic-parabolic regularization, con¬ 
dition V (along with (H0)-(H3) and (Al)-(A3) below) has been shown suf- 
hcient for establishing nonlinear stability for Lax and overcompressive shock 
prohles |Ral IMaZ.41 IZ.4j . Except in the case of |Raj , these nonlinear analyses 
are carried out through consideration of the perturbation 

v{x, t) = u{x, t) — 

where 6 {t) is introduced as a local phase chosen to insure that the shapes of u 
and u are compared, rather than their positions (in this way, orbital stability 
is considered). In this context, and under the assumption of stability criterion 
{T>), it can be shown that for initial perturbations 

|n(a:, 0)1 < Eo(l \x\)~^^‘^, 
some Eq sufficiently small, there holds 
(1.6) \\v{x,t)\\LP < CEo{l + t)~^^^~p\ 

from which we observe asymptotic decay in time for all p > 1 . 
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A natural refinement of this type of analysis regards the consideration 
of diffusion waves, which are defined as exact solutions to a family of con- 
vecting Burgers’ equations (convection along outgoing characteristics of the 
underlying hyperbolic problem), and which carry precisely the mass in 
(USD that does not decay asymptotically in time |Lin21 ILindj . Considering, 
then, the perturbation 

v{x, i) = u{x, t) — u^^^\x) — ip{x, t), 


where (p{t,x) represents the sum of diffusion waves, we have 


f + CXD 


b + OO 


if>{x,ff)dx= I {u{x,0) — u^^^\x))dx, 


and consequently. 


^ + 00 


v{x, 0)dx 


0 . 


In this way, we reduce to the case of zero-mass initial data, for which pertur¬ 
bations decay equally from above and below the profile, and the asymptotic 
rate of decay is doubled (at least at the linear level, and for an initial pertur¬ 
bation with sufficiently fast spatial decay). Working in this context, Raoofi 
has shown that stability criterion {T>), along with (H0)-(H3) and (Al)-(A3), 
are sufficient for establishing that for initial perturbations 


v{x, t) 


u{x, t) — u^* {x) — ip{x, 0) 


du^ 


m, 


where 5* is the asymptotic shape/location of the shock, and with \u{-, 0) —u 
sufficiently small, there holds jRaj 

|k(-,t)||LP<CEo(l + t)-5(l-F)-i 


The goal of the current analysis is both to refine the analysis of Raoofi to 
the case of pointwise (rather than L^) estimates, and to introduce a simplified 
bootstrapping argument through which estimates on the perturbation v{t, x) 
emerge in straightforward fashion. Our estimates are similar to those of Liu, 
developed for Lax type prohles in the case of identity viscosity and under the 
assumption of weak shock strength. Our analysis has no such limitations, 
though we note that the form of our viscosity and the shock strength under 
consideration are encoded in our spectral assumptions. 
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Throughout the analysis, we will work in a coordinate system moving 
along with the shock, so that without loss of generality, we consider a standing 
prohle u{x), which satishes the traveling-wave ordinary differential equation 
(ODE) 

(1.7) B{u)u = F{u) - F{u_). 

Considering the block structure of B, this can be written as: 

( 1 . 8 ) F\u^, = F\ui, 

and 

(1.9) bi{u^y + b2iu^^y = F^\uyu^^) - F^\uy,u^J). 

We are interested in the asymptotic and pointwise behavior of ft, a solution 
of p. 2 |l and a perturbation of ft. 

We assume that, by some invertible change of coordinates u w{u), 
possibly but not necessarily connected with a global convex entropy, followed 
if necessary by multiplication on the left by a nonsingular matrix function 
S{w), equations (ll. 2 |l may be written in the quasilinear, partially symmetric 
hyperbolic-parabolic form 

(1.10) A^wt +Awx = {Bwx)x +G, w = 

G G W, X G M'^, f G M, where, dehning W± := W{U±): 

(Al) A{w±), All, A^ are symmetric, A^ > 0 . 

(A2) Dissipativity: no eigenvector of dF{u±) lies in the kernel of B{u±). 
(Equivalently, no eigenvector of A(A°)“^(tt;±) lies in the kernel of B{A^)~^{w±).) 
(A3) 5 = , G = , with Reb{w) > 9 for some 9 > 0, for all W, 

and g{wx,Wx) = OdtCxP)- Here, the coefficients of dl.lOj] may be expressed 

in terms of the original equation (O, the coordinate change u —>• w{u), and 
the approximate symmetrizer S{w), as 

A° := S{w){du/dw), A := S{w)dF{u{w)){du/dw), 

^ ^ B := S{w)B{u{w)){du/dw), G = —{dSwx)B{u{w)){du/dw)wx- 
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Along with the above structural assumptions, we make the technical hy¬ 
potheses: 

(HO) F, B,w, S e C®. 

(HI) The eigenvalues of A* := Aii(A'j’]^)“^ are (i) distinct from 0 ; (ii) of 

common sign; and (iii) of constant multiplicity with respect to u. 

(H2) The eigenvalues of dF{u±) are real, distinct, and nonzero. 

(H3) Nearby u, the set of all solutions of (ll.lj) - (ll. 2 |l connecting the same 
values u± forms a smooth manifold {h'^}, 6 eU G = it. 

We note that structural assumptions (Al)-(A3) and technical hypothe¬ 
ses (H0)-(H2) are broad enough to include such cases as the compressible 
Navier-Stokes equations, the equations of magnetohydrodymics, and Slem- 
rod’s model for van der Waal gas dynamics jZ.4j . Moreover, existence of 
waves u satisfying (H3) has been established in each of these cases. 


Definition 1.1. An ideal shock 


( 1 . 12 ) 


u{x, t) 


U- X < st, 
M+ X > st, 


is classified as undercompressive, Lax, or overcompressive type according as 
i — n is less than, equal to, or greater than 1, where i, denoting the sum of 
the dimensions and of the center-unstable subspace of df{u-) and the 
center-stable subspace of df{u+), represents the total number of characteris¬ 
tics incoming to the shock. 

A viscous profile (HU is classified as pure undercompressive type if the 
associated ideal shock is undercompressive and i = 1, pure Lax type if the 
corresponding ideal shock is Lax type and i = i — n, and pure overcompressive 
type if the corresponding ideal shock is overcompressive and i = i — n, i as 
in (H3). Otherwise it is classified as mixed under-overcompressive type; see 

\lM\z^. 

Pure Lax type prohles are the most common type, and the only type 
arising in standard gas dynamics, while pure over- and undercompressive 
type prohles arise in magnetohydrodynamics (MHD) and phase-transitional 
models. 

Under assumptions (A0)-(A3) and (H0)-(H3), or their analogs in the 
real viscosity case, condition ifD) is equivalent to (i) strong spectral stability, 
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ct(L) C {3?A < 0}U{0}, (ii) hyperbolic stability of the associated ideal shock, 
and (iii) transversality of h as a solution of the connection problem in the asso¬ 
ciated traveling-wave ODE, where hyperbolic stability is dehned for Lax and 
undercompressive shocks by the Lopatinski condition of |M.llIM.21 IM.31 IFrej 
and for overcompressive shocks by an analogous long-wave stability condition 
(see iVii) below); see [ZHl IMaZ.31 IZSl ITT^ IZ.31 IZ.4j for further explanation. 
From now on, we assume {T>) to hold along with (A0)-(A3) and (H0)-(H4). 
We also assume that the the shock is a pure Lax or overcompressive one. 

Setting A± := df{u±), r± ;= cPf{u±), and B± := B{u±), denote by 

(1.13) < 02 <■■■< a~ and af < < ■ ■ ■ < 

the eigenvalues of A_ and A+, and if, rf left and right eigenvectors as¬ 
sociated with each af, normalized so that {ljrk)± = Si, where 51 is the 
Kronecker delta function, returning 1 for j = k and 0 for j ^ k. Under this 
notation, hyperbolic stability of u, a Lax or overcompressive shock prohle, is 
the condition: 

(Dii) The set {r=‘=; ^ 0} U ^ = 1,''' 5 forms a basis for 

M"', with computed at ^ = 0. 

As said before, (Dii) is satished whenever (D) holds. 

Dehne scalar diffusion coefficients 

(1.14) lif := 

and scalar coupling coefficients 

(1.15) Ty=(!jr(r„r;))±. 

Following jLiull ILin21 ILin3j . dehne for a given mass rrij the scalar dif¬ 
fusion waves ip~{x,t-,m~) as (self-similar) solutions of the Burgers equations 


( 1 - 16 ) = -%■ )')- 

with point-source initial data 

(1.17) (p“(x, -1) = mj5o{x), 
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and similarly for Given a collection of masses prescribed 

on ontgoing characteristic modes aj < 0 and > 0, define 

(1.18) (p{x,t) = ipj{x,t]mj)rj + Y 

a~ <0 a+ >0 

Also define 


(1.19) 


and 


( 1 . 20 ) 


'ipi{x,t) := x{x,t) Yi^^^‘^{^ + \x - ajt\+t3) 

a-<o 

+ X{x,t) ^(l + f)“^/^(l + |x-a+f| +^^)■^/^ 
a +>0 


'i/ji{x,t) := x{x,t) Yi^ + ^) ^-^^(1 + 1^: - a^Gl) 

a-<o 

+ X{x,t) ^(l + f)“^/2(l + |a;-a+^|)■^/^ 

a +>0 


where x{x, f) = 1 for x G [a^ t, a'^t] and zero otherwise. Also, 

^p2{x, t) := Y^ (1 + |x - a~t\ + 

a ~<0 

( 1 . 21 ) 

+ 5 ](l + |x-a+t|+t^/2)-3/2^ 

a +>0 

and 

( 1 . 22 ) a(x, t) := x(x, f)(l + t)“^/^(l + |x|)“^/^ 

Theorem ll. 2 l and Corollaries 1 1. 41 and 1 1. 51 are the main resnlts of this paper: 

Theorem 1.2. Assume (Al)-(A3), (H0)-(H3) and iV) hold, and u is a 
pure Lax or overcompressive shock profile. Assume also that u solves 03) 
with initial data Uq and that, for initial perturbation Uq := Uq — u, we have 
\uo\L^nm < Eq, |Mo(a^)| < ^o(l + \x\)~^, and \d,,Uo{x)\ < Eo{l + |x|)“ 2 , for 


Eq sufficiently small. Then there are an i-array function 6(t), and a small 
constant 5^, such that if v := u — u * — (p — ^5, then 

(1.23) |u(a;, t)| < CEo('0i +'02 + a), 
and 

(1.24) \vffx,t)\ < CEot~^{l + t)5(^/!i + 02 + a); 
furthermore, 

(1.25) \6{t)\<CEo{l + t)-"2, 
and 

(1.26) \5{t)\<CEffl + t)-\ 

for some constant C (independent of x,t and Eq). 

The proof of Theorem 1 1.21 uses a straightforward bootstrapping argument, 
combined with estimates proved in jR,aj and restated here in the following 
proposition. 

Proposition 1.3. Under the conditions of Theorem Al.iA 

(1.27) H;t)\LP < C'Eo(l + f)-5(i-F)-i, 
for any p, I < p < oo; also 

(1.28) \v{-,t)\Hs < CEq{ 1 + t) 2^ 
for any s < 4, and for some constant C. 

Proof of Proposition II.3L See jRaj. In |Raj the bound ()1.28j) was es¬ 
tablished for s < 3; basically, to prove (HHHI), it was shown, using energy 
estimates, that the i7^-norm of u is controlled by the i7^-norm of the initial 
data and the L^-norm of u. However, the same proof (with a slight modifi¬ 
cation) can be used for s = 4; we need only to assume that the initial data is 
small in and the coefficients are in (rather than as in jR aj ). The 
reason we need two times differentiability of the coefficients is in the fact that 
we need the exponential decay of u and its derivatives to their endstates up 
to 2s derivatives; see (5.56) in |Raj and Lemma f2. II in the present work. □ 
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Taylor expansion gives ns 


(1.29) + 0(|5(t)| 

do 

since it approaches in an exponential rate to its endstates (see Lemma in}. 

Bnt is then smaller than the right hand side of (ll.2dll . Hence we 

have the following. 

Corollary 1.4. Under the assumption of Theorem M/A 

(1.30) u{x, t) — = ip{x, t) + + "02 + tt). 

Also, 

(1.31) \u{-,t) - - ip{-,t)\LP = 0((1 + f)“2d-p)-4). 


Withont “instantaneous shock tracking” d{t), however, we obtain the 
following, which Liu proved for the artificial viscosity case |Lin3j . 

Corollary 1.5. Under the assumption of Theorem M.A 


(1.32) \u{-,t) - U^* - ip{-,t)\Lv 


0 ((l + f) 2 d p) 4 ) forl<p <2 
0((l + t)“5) for2<p<cx). 


The picture of asymptotic behavior described in Theorem II .21 and corol¬ 
lary II .41 was introduced on heuristic grounds by Liu [Tln2j in the context 
of small-amplitude Lax-type shock waves and artihcial (identity) viscosity 
B = I, and, along with the accompanying analysis of |Linlj described be¬ 
low, played an important role in the subsequent analysis by Szepessy and 
Xin in jSXj establishing for the hrst time stability (with no rate) of small- 
amplitude Lax-type shock prohles with B = I. Bounds i fTTri ) validating this 
picture were established for large-amplitude Lax or overcompressive prohles 
and general viscosity by Raooh [Raj. See also earlier pointwise arguments 

Liu3[ 


sketched (but not completed) 


m 


In the proof of Theorem 11.21 


we will assume and use the bounds (071) and (USB}, already established in 
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Remark 1.6. The estimates of Theorem o are similar to, though weaker 
than, those developed by Liu in the case of identity viscosity and weak shock 
strength jLindj . In particular, Liu’s estimates provide slightly sharper results 
in terms of the pointwise bounds, and also give more information along the 
characteristic modes. In order to provide results as sharp as those of Liu, 
we would have to treat the critical nonlinearity (pv similarly as we treat the 
nonlinearity 93 ^, which we analyze by a clever approach of Liu’s in which 
he integrates the nonlinear interaction integral by parts in t, making use of 
the observation that {dt + a~dx)(pj decays similarly as the t-derivative of 
a heat kernel |Linlj . (For the argument in our context, see ()4.d7j) and the 
surrounding discussion, for which under a shift of coordinates, differentiation 
with respect to r replaces the operator {dt + ajd^)- The resulting estimate is 
stated in of Lemma fd.21 1 In order to apply this approach to the term 

ipv, we would additionally have to carry an estimate on {dt + a~dx)v through 
our argument. We leave the full details of this calculation to a separate 
paper [HRZ] . We remark also that the estimates of Liu are uniform in shock 
strength e := |m+ — M-| as e —>• 0. More specihcally, e appears explicitly in 
Liu’s estimates, and the estimates remain largely unchanged as e ^ 0, though 
one of the estimates (denoted Xi) increases slightly due to the loss of a term 
as a possibility in an estimate that takes the minimum of three 
quantities (the other two quantities remain uniform in e) jLiu3j . In our case, 
shock strength is assumed fixed, and our estimates are not uniform as e —0. 
In particular, coefficients in the Green’s function estimates we use blow up 
as e —^ 0, and we must counter this by reducing Eq, the small constant 
multiplying our initial perturbation. We regard the extension of our analysis 
to the case uniform in (small) shock strength as an important future project. 
One of the advantages of our approach, on the other hand, is the use of the 
instantaneous shock tracking 6{t), as well as using the previously established 
norms, which makes our proof both more straightforward and easier to 
generalize. We also remark that since we proceed from Green’s function 
estimates obtained for possibly large-amplitude shock prohles (information 
about the amplitude is encoded in the spectrum of the linearized operator), 
our analysis applies to this case as well. Finally, the estimates of Gorollary 
o recover the estimates of |R aj . with a slight improvement on the shift 
location estimates \S{t)\ and S{t). As mentioned in jR aj . these estimates, in 
the case p > 2, are a slight improvement of Liu’s. 


Remark 1.7 (Remarks on the Strictly Parabolic case). The case 
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of the strictly parabolic systems can be treated in a very similar way, the 
proof being almost identical. However, we need fewer assnmptions for the 
eqnations, or for the initial pertnrbation, in this case. Basically, instead of 
(A1)-(A3) and (H0)-(H3), we assnme the following assnmptions in the case 
of strictly parabolic case. 

(HO) f,BeC\ 

(HI) Re cr{B) > 0. 

(H2) The eigenvalnes of df{u±) are real, distinct, and nonzero. 

(H?>) For some 0 > 0, and all real /c, we have 

Rea{—ikdf{u±) — k‘^B{u±)) < —Ok'^. 

(H4:) The set of all stationary solntions of ()l.l|l - ()1.2j) near it, connecting 
the same valnes u±, forms a smooth manifold {h^}, 6 G M^, it^ = it. 

We also assnme that (D) holds. For the initial perturbation we only need 
Uq G for some 0 < a < 1. In Theorem 11.21 concerninEf the real viscosity 
(partially parabolic) case, we require the pointwise bound on the derivative 
of the initial perturbation—i.e. the bound |c?a;no(3^)| < -Fo(l + —only 

in order to control the derivative in the hyperbolic modes, which are absent 
in the strictly parabolic case. The necessary bounds on the derivatives in the 
strictly parabolic case can be achieved through short time estimates similar 
to ones carried out in On the other hand, as we are assuming 

less on the initial data, does not necessarily hold. See |R,aj for more 

details. 


2 Linearized equations and Green function 
bounds 

Before stating the Green function bounds for the linearized equation, we 
need some preparation to do. First we need the exponential decay of it to its 
endstates. The following lemma proved in jMaZ.2j provides us with that. 

Lemma 2.1. Given (H1)-(H3), the endstates u± are hyperbolic rest points of 
the ODE determined by (US) on the r-dimensional manifold i.e., the 
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coefficients of the linearized equations about u±, written in local coordinates, 
have no center subspace. In particular, under regularity (HO), 

(2.1) DiDi{u\x)-u^) = 0{e-‘^\^\), a > 0, 0 < j < 8, z = 0,1, 
as X —>■ ±cx). 

Instead of linearizing abont m(-), we linearize abont where 5*, de¬ 

termined a priori by the mass of the pertnrbation, would be the asymptotic 
location or shape of the shock. Linearizing around gives us 

(2.2) vt = Lv := -{Av)x + {Bv,c)x, 
with 

(2.3) B{x) ;= B{u^*{x)), A{x)v := df{u^*{x))v — dB{u^*{x))vu)f. 
and let 


(2.4) G{x,t-,y) := e^^6y{x). 

Denoting A^ := y4(±cxo), B^ := B{±oo), and considering lemma ITTl it 
follows that 


(2.5) \A{x) - A-\ = 0{e-^\^), |5(x) - 5”| = 0(e-^l"l) 


as a: —>■ —oo, for some positive r]. Similarly for A~^ and B~^, as x —>■ -|-oo. Also 
|A(x) — A^\ and \B{x) — B^\ are bounded for all x. 

Define the (scalar) characteristic speeds af < ■ ■ ■ < af (as above) to be 
the eigenvalues of A^, and the left and right (scalar) characteristic modes if, 
rf to be corresponding left and right eigenvectors, respectively (i.e., A^rf = 
afrf, etc.), normalized so that If ■ rf = 61 and If ■ r^ = 5^. Following 
Kawashima |Kawj . define associated effective scalar diffusion rates (lf\j = 
1 , • • ■ , n by relation 


( 2 . 6 ) 


(Pt 0 \ 

vO 


= diag L^B^R^, 


where := {if,..., iff, := {rf ,..., rf) diagonalize A^. 
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Assume for A and B the block structures: 


A 


\^21 ^22 / 




Also, let a*(a;), j = 1,..., (n — r) denote the eigenvalues of 

A* := An — A12B22 B21, 


with l*{x), r*{x) e M” ^ associated left and right eigenvectors, normalized so 
that = 61- More generally, for an m*-fold eigenvalue, we choose {n — r)x 
m* blocks L* and R* of eigenvectors satisfying the dynamical normalization 

Lfd,R* = 0, 

along with the usual static normalization LfRj = dllm*', as shown in Lemma 
4.9, |MaZ.lj . this may always be achieved with bounded L*, R*. Associated 
with L*, R*, dehne extended, n x m* blocks 





Eigenvalues a* and eigenmodes £*, TZ* correspond, respectively, to short-time 
hyperbolic characteristic speeds and modes of propagation for the reduced, 
hyperbolic part of degenerate system m- 

Dehne local, mj x mj dissipation coefficients 


V*{x) :=-L* D^R*{x), j = 1,..., J < n 


where 


D*(x) := 


A12B22 



A 22 B 22 B 21 + A*i?22^-B21 + B22dx{B22 B 21 ) 


is an effective dissipation analogous to the effective diffusion predicted by 
formal, Chapman-Enskog expansion in the (dual) relaxation case. 

At a; = ±cx), these reduce to the corresponding quantities identihed by 
Zeng IZilllLZi] in her study by Fourier transform techniques of decay to 
constant solutions U = u± of hyperbolic-parabolic systems, i.e., of limiting 
equations 

Ut = L±U := —A±Ux + B±Uxx- 
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As a consequence of dissipativity, (A2), we obtain (see, e.g., jKawL ILZel 
IMaZ.3j i 


(2.7) f3f > 0, Rea{ri*^) > 0 for all j. 

However, note that the dynamical dissipation coefficient D^{x) does not agree 
with its static counterpart, possessing an additional term i?225x(-B72^i?2i), 
and so we cannot conclude that (EID holds everywhere along the prohle, 
but only at the endpoints. This is an important difference in the variable- 
coefficient case; see Remarks 1.11-1.12 of jMaZ.3j for further discussion. 

Proposition 2.2. IMaZ. Under assumptions (Al)-(A3), (H0)-(H3), and 
{V), the Green distribution G{x,t]y) associated with the linearized evolution 
equations may be decomposed as 

(2.8) G{xG]y) = H + E + S + R, 
where, for y < 0.- 


(2.9) 




J 

■=^a*j~^{x)a*{y)n*{x)C*{y,t)6a:-a*t{-y)Rf{y) 

1=1 

J 

1 = 1 


where the averaged convection rates d* = d*{x,t) in 113. .91) denote the time- 
averages over [0, t] of a*{x) along backward characteristic paths z* = z*{x,t) 
defined by 

dz*/dt = a*{z*), z*{t)=x, 

and the dissipation matrix C,* = (*{x,t) G is defined by the dissipative 

flow 

dQ/dt = -y*{z*)Q, C;(0)=/™,; 

and Sx-ap denotes Dirac distribution centered at x — d*t. 


( 2 . 10 ) 


E(x,t-, y) 


I 


E 


du^{x) 


ei(2/,i), 

15 = 5 * 
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(2.11) ej{y,t) ;= ^ I errfn (] — errfn i ^ | | / 


a, >0 


V V^l^k ^ 


V^f^k ^ 


jfc 


( 2 . 12 ) 

^(a;,t;|/) := X{t>i} Y1 C'(47r/3fc ' 


aj. <0 


“fc >0 


gX g-X 


+ X{t>i} Y1 ^k\^^Pjkt) 

“fc >0, a“ <0 

a^>0, af>0 


\y\ 


gX g-X 


gX _|_ g X 


with 


(2.13) 

4k(y^ ■= ( 

and 


(2,14) 

■= + 


bl /«, 



(^k 1 

The remainder R and its derivatives have the following bounds. 

(2.15) 

fi(i,«;t/) = 0(e-’<l‘'-»l+‘>) 

n 

+ ^ O ((t + l)-^/‘^e-^^^ + ^-lt2^-{x-y-altflMt 

k=\ 

>0, a~ <0 

+ E ^11% i|>l.llO((l + l)-'/2t-i/'2)e-(—;«-l»/<>:l))V*«e-i-- 


aj, >0, a+ >0 
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(2.16) 

./ 

R,(x, t,y) = Yl 0{e-'‘)S,-,;,{-y) + 0(e-’><l"’-"l+6) 

i=i 


n 

+ ^ O ((t + l)-^/2g-r7X+ ^-ri\x\ /Mt 


k=l 




a^, >0, a^. <0 




—■qx 


dk >0, at >0 


(2.17) 

J 

Rx{x,t-,y) = J20{e-^^)6x--ap{-y) + 0(e-’'(l"-^l+*)) 
i=i 


n 

+ ^ O ((t + t-i(t + i)i/2g-(^-i/-a^t)- 


/Mt 


fc=l 




—qx^ 


dk > 0 , dj <0 




—qx 


dk >0, at >0 


Moreover, for \x — y\/t sufficiently large, |G| < Ce as in the 

strictly parabolic case. 

Setting G := S + R, so that G = H + E + G, we have the following 
alternative bounds for G. 


Proposition 2.3 ( [ZHJ . |MaZ.3i IHZJ I. 


Under the assumptions of Propo- 
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sition\2^ G has the following bounds. 

\d^,yG{x,t-,y)\ < 

n 


k=l 


(2.18) 


I V h-\ylH \))'^ 

+ ^{Kt\>\y\}i e J e , 


“fe >0>aj <0 


I V +-l/2„-(ai-a+(4-|y/afc |))^/Mt -r;x ^ 


>0, aj >0 


for y < 0, and symmetrically for y > 0, for some rj, C, M > 0, where af are 
as in Provosition \2.iA (3^ > 0, denotes the positive/negative part of x, and 
indicator function X{\a-t\>\yW 1 \^kA — Ivl 0 otherwise. Moreover, 

all estimates are uniform in the supressed parameter 5*. 


Remark 2.4. We will refer to the three differently scaled diffusion kernels 
in (ITTSl) respectively as the convection kernel, the reflection kernel, and the 
transmission kernel. We recall the notation 


errfn( 2 ;) 



e ^^df. 


Remark 2.5. Green function bounds for the strictly parabolic case (see 
Remark EZD is very similar. The main difference is that, in the strictly 
parabolic case, the hyperbolic part H is absent in the decomposition of G as 
in (ITHl) . 

Remark 2.6. From we obtain by straightforward calculation (see 

jMaZ.dj i the bounds 


(2.19) 


\dtej{y,t)\ < Gr^!’^ ^ ^-\y+ait\pMt^ 

a->0 

\dyej{y,t)\ < Gr^/^ ^ ^-\y+a^t\GMt 
a^>0 

\dyte,{y,t)\ < Gt-^ 


errfn 


y-^kt 

VWt 
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for 1 / < 0, and symmetrically for ?/ > 0. 


3 Nonlinear analysis 

Let u solve (O, and, using ("Pii), assume that 



- u{x) = 


mjVj 


■<o 


E 

a +>0 


rrijrj 



|<s=o 


with rriiS and Cj’s small enough. Using the Implicit Function Theorem, we 
can hnd 5* such that 



u^* 


.X = 




mf. 


a. <0 


aj>0 


where each m' is just “slightly” different from rrii. Notice that this way 
we have no “mass” in any / ^ direction anymore. With a slight abuse of 
notation we drop the ' sign from each m' and denote it simply by m*. 

Remark 3.1. In the case of Lax-type shock waves, u^{x) = u[x + h), hence 
^ = u'{x)^ and h* can be explicitly computed: h* = ci. 

Let u{x,t) = u{x,t) —u^*[x) and use Taylor’s expansion around u^*{x) to 
hnd 


(3.1) ut + {A{x)u)^ - {B{x)ux)x = -(r(a;)(M, u))x + Q{u, Ux)x, 
where T{x){u,u) = cPf{u^*){u,u) — (PB{u^*){u,u)ul* and 

Q{u,Ux) = 0(|M||Ma,| + \uf). 

Denote T^*" = r(±cx)). Dehne constant coefficients bfj and Tj^-^ to satisfy 

n n 

(3.2) = = 

i=l i=l 

Then, of course, (5f = bf^ and yf := Tf--. 
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Now define by (jl.l(i|l and (ITTTIl . and likewise for ipf. Finally ip is 
defined by (jl.lSj) . Then set v := u — ^p — with 6{t) to be dehned 

later, and assuming 5(0) = 0. Notice that, by our choice of 5* and diffusion 
waves ip^^s, we have zero initial mass of v, i.e.. 


(3.3) 


^+oo 


v{x, 0)dx = 0 


Replacing u with n + 9 ? + ^S{t) in ()3.1|) ( ^ computed at 5 = 5*), and 

using the fact that ^ satishes the linear time independent equation Lv = 0, 
we will have 


(3.4) 

where 

(3.6) 


du^ du^ • 

Vt- Lv = <F(x, t) + J^((p, V, 


= o(|.„|^ + bllH + |,. 11^61 + |„||^6| + 1 ^ 61 '^ 


86 


86 


86 


86 

+ \{(p + v + -—6){(p + V + wtw<5)x| + \p> + v + 5|^). 


86 


86 


86 


Furthermore 


f)ll^ 

J^{v,^,-^6{t))^ = 0{T{v,(p,—6) 


(3.6) 


+ \i^’ + L + ^S)^\\{v^^ + (p + ^6)^\ 

+ \'>^ + L + + p} + ^6)^^\), 


86 ^ 86 

and <F(a:, t) := -p>t - {A{x)p>)x + {B{x)p>x)x - (r(a:)((^, p)))x- For <F we write 


(3.7) 

<h(a;, t)=- ((pi + A(p:^ - Bp>^^ + F((p, pi)^) 

= - Y1 T'X" + iMx)p^'r-)x - iB{x)p>lr;)x + (F((r)((pVr, (^V,"))^ 

a- <0 

- Y1 - {B{x)(plrt)x + iT{x){(pWl,(p^rl))x 

a+>0 

- '^iLiP>jPix){r^,rf))x- 
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Let us write a typical term of the first summation (a^ < 0)in the following 
form: 

^3-8) 

- {B{x)^lr-)^ + (r(a;)((pV-, (pV”))^ 

= - {B{x) - B-)(plr; + (r(a:) - T-){(pW-, cpW;)]^ 

+ - i^l^B-r-) + 

Now we use the dehnition of (p* in and the dehnition of coefficients bij 

and Tijk in (jd.2j) to write the last part of (jd.Sjl in the following form: 


vVi + r- )-i^i^B r- ) + (r. ,ri )) 

=-v’LE'-.WE r 

j¥=i 


r ■ 

3 


Similar statements hold for > 0 with minus signs replaced with plus signs. 
Now we employ Duhamel’s principle to get from itra 


(3.10) 




v{x,t) = / G{x,t-,y)vo{y)dy 


' —OQ 

f*t /• + 00 



du^ 


+ I I G{x,t - s-,y)B{ip,v,—6)y{y,s)dyds, 


' 0 J —oo 

(*t /•+O0 



G{x, t — s] y)^{y, s)dyds 


'0 J — oo 

du^ 


To obtain the above, we used 

r+oo 


,du\ . , Lt^U^ 

G(.,t;y) — (y)dy = e jr = JT 


and (5(0) = 0. 
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Set 


m = 


(ii{y,t)vo{y)dy 


— OO 

t /• + 00 


(3.11) 



du^ 


' 0 J —CO 
r+OQ 



ei{y, t - s)J^{ip, V, —6)y{y, s)dyds 
ei{x,t- s;y)^{y,s)dyds. 


' 0 J —CO 

Using iHI) and G = H + E + G we obtain: 

/ +00 

{H + G){x,t-,y)vo{y)dy 

■OO 

{H + G){x,t - s-,y)E{^,v,—6)y{y,s)dyds 

pt /‘H-oo 

/ / {H+ G){x,t - s-,y)^{y,s)dyds. 

'o J -OO 


(3.12) 


+ 

+ 


In addition to v{x,t) and d{t), we will keep track in our argument of 
Vx{x,t) and 6(t), the latter of which satisfies 


m = 


dtei{y,t)vo{y)dy 


— OO 

t r+oo 


(3.13) 



du^ 


' 0 J —OO 
/• + 00 



dtei{y, t - s)E{ip, v, s)dyds 

dtei{y,t- s)^{y,s)dyds, 


0 J —CO 


where we have taken advantage of the observation, apparent from 
that e{y,0) =0. 

The following lemmas are the main ingredients for the proof of Theorem 
11.21 Their proofs, however, are postponed to Section HI 

Lemma 3.2 (Estimates for linear part). Suppose that for some Eq > 0, 
we have that Vo{y) satisfies the conditions 

\v„(y)\ < Eo(l + \y\)-^'^ 

r+oc 

/ Vo{y)dy = 0. 
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Then there holds 

(3.14) 

(3.15) 

(3.16) 


^+oo 


G{x, t] y)vo{y)dy < CEo'ijj 2 {x, t). 


'‘+O0 


J —OQ 
'H-oo 


ei{y,t)vo{y)dy <CEo{l + t) 
dtCiiy, t)vo{y)dy < CEo{l + 


where C does not depend on Eq. 

3 

Lemma 3.3. If |fo((z:)| < Eo{^ + 1^1)”^) then 


r‘+oo 


(3.17) / H{x,t-,y)vo{y)dy < CEoe + \x\) 2 

J —00 

If\d^Vo{x)\ < Eo{l + \x\)~T then 

p +00 ^ 

(3.18) / H^{x, t; y)vo{y)dy < CEoe~^\l + |x |)“2 


for some 9 > 0. The right hand sides of and are obviously of 

order a{x, t). 

Lemma 3.4 (Estimates for nonlinear part). For G{x,t-, y) as in Propo- 
sitionsW~^ a,nd,Wf^ we have 


t /• + 00 


(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 



\Gy{x,t- s;y)\^{y,s)dyds < G{iJi + -02 + a)(x,t), 


0 J —00 
t /•+00 



G{x,t- s;y)^{y,s)dyds < GEo{f)i + ii 2 ){x,t), 


' 0 J —00 
f*t /•H-oo 



\dyei{y,t- s)\'^{y,s)dyds < G{l + t) 


0 J —00 

t /•+00 



ei{y,t - s)^{y,s)dyds <GEo{l + t) 


0 J —00 
t /•+O0 



dtei{y,t- s)^{y,s)dyds <GEo{l + t) \ 


0 J —00 
t /*+c>o 



\dytei{y,t - s)\'I/{y, s)dyds < G{l + t) \ 


0 J—00 
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for $(?/, s) as in (O and 


(3.25) 


'^{y,s) = {1 + s) ^^‘^{f)i + 'ijj2 + a + ip){y,s) 


Lemma 3.5. // |T(?/, s)| < s + 1^2 + a + (p){y, s) + s then 

f*t /*+O0 


(3.26) 



H{x, t — s] y)T{y, s)dyds < C{'tpi + V ’2 + a){x, t). 


'0 J —00 

If \dyT{y, s)| < s“^/^('^i + 1^2 + <y + (p + e~^^y^){y, s), then 

f*t /•+00 


(3.27) 
Also, 

(3.28) 



H,,{x,t - s-,y)T{y,s)dyds < C(V’i + ^2 + a)(a;,t). 


0 ^ —00 


t /•+00 



H{x,t — s; y)^{y, s)dyds < CE^i^i + V ’2 + a){x, t){x, t), 


' 0 j —OQ 

and similarly with H replaced by 


Proof of Theorem EH We use estimates and in order to 

obtain a supremum norm on v and its first, second and third derivatives in 
X. We prove (II2SD with -01 replaced by -^iji.e. we prove 


(3.29) \v{x,t)\ <CEo{'ljJl+^jJ 2 +a);, 


3 

but then, as we have |n(a;, t)| < CEo{l + t)~i (this is p.27jl for p 
as 

3 — 

min{(l + t)“4, 'ijj^{x,t)} ~ i>i{x,t), 


00 ), and 


the proof of would be immediate. 

Let 


(3.30) 

m ■■= 


sup 

y,0<s<t 


\v{y,s)\ ^ _ |w(y,s)| _ 

(V’l + 1/^2 + a) ( 2 /, s) y,o<s<t t“i(l + t)^V’i + -02 + Ce){y, s) 
+ sup |(5(s)|(l + s)5 + sup |(5(s)|(l + s). 

0<s<t 0<s<t 
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Our aim is to show that 


(3.31) at)<C{Eo + Eoat)), 

from which we conclude ({t) < if Equivalent to (Id.dlj) 

is 

(3.32) \vix,t)\ < C{Eo + EoC{t)){^i + i >2 + a){y, s), 


(3.33) m\<c{Eo+Eoami+t)-K 

and similar statements for and <5. 

Estimates ofv{x,t). Looking at (I3.12j] . there are three parts that should be 
estimated. The hrst part, the linear part, is carried out using (IXT7I) and 
dsm. For the second part, we notice that, by (HH), (ESDI), OH), (11123), 
(11123, Lemma I2II1 and the dehnition and bounds of </9, 


|.i^(l/,s)| < C{Eo + C{t))^{y,s) 

with E{y, s) as in (I3.5|l . and T as in ()3.25|1 . Hence, by integration by part. 


/■* ~ du^ 

/ / G{x,t - s]y)E{(p,v,—6)y{y,s)dyds 

G ~ du^ 

/ / Gy{x,t-s]y)E{^p,v,—6){y,s)dyds 

< G{Eq + EoC{t))h{x,t), 


by (13.191) . 

As we do not have good estimates for Hy we cannot do the same with 
part containing H, so, instead, we notice that 

\Ey{y,s)\<G{Eo + at))Tiy,s) 

Ey{y, s) as in (I3.b|l . and T{y, s) as in Lemma ITHl (The subtle fact is that, as 
we have Vxx in tFy{y^ s), we need to use (ll.28|l . hence having only (1 + 1)~^, 
instead of (1 + f)“4 which appears in T.) Therefore we can use (j3.2(H) to 
obtain the desired estimate. The third integral of (IXT^ can be estimated 
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similarly using innnii and (HOHll . 

Estimates of 6(t) and 6{t). Similarly and using (Id.111) and (ld.21Kld.24|l . 
Estimates ofVx{x,t). By (Id.l2|l . 


^H-OO 


Vx{x,t) = 

(3.34) + 

+ 


{H^ + Ga:){x,t;y)vo{y)dy 


’ —OO 

f*t /• + O0 



du^ 


' 0 J —OO 
f*t /• + 00 



{H^ + Ga:){x, t - s; y)J^{ip, v, —6)y{y, s)dyds, 
+ G^){x,t- s;y)^{y,s)dyds. 


0 J —OO 


The parts involving are treated using (ITTHl) . (IT?rl) and the similar state¬ 
ment for <h in Lemma ItTHl this because, due to ()1.28D Tyy{y, s) has the bounds 
similar to iFy{y, s). For Gx {x,t;y)vo{y)dy, we notice that Gx is at least 
as good as G away from 0. Therefore, for f > 1, a proof identical to that of 
(ITTTl) can also prove 


“+CXD 


Gx{x,t;y)vQ{y)dy < GEo'ip 2 {x, t). 


On the other hand, for f < 1 we use the fact that Gx{x, t; y r\j t 2G{x,mt-,y) 
for t near zero and some positive m. This, together with (It:^ provides us 
with the necessary bounds for the linear part of the calculations. 

F inally. 


* /■+“ ~ du^ 

Gx{x,t- s]y)E'{ip,v, —6)y{y,s)dyds 



0 J —OO 


+ 


85 

r*-i r+^ _ 

I Gxy{x,t - s-,y)J^{ip,v,-^5){y,s)dyds 

8vf 

Gx{x,t- s]y)J='{ip,v, —5)y{y,s)dyds. 


J 0 J —OO 
pt p-\-oo 

h-1 J-OQ 


The hrst integral can be computed exactly the same way (imi is proved, 
for away from t — s = 0, Gxy{x, t — s;y) is as well as Gxy{x, t — s;y). For the 
second integral, a close investigation of the proof of (ixn shows us that, if 
we limit the time integration to t — 1 < s < f, then the same proof works to 
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prove that 



s;y)\\T{y,s)\dyds < C{ijji + ^|J 2 + a){x,t), 


with T as in Lemma 13.51 This fact plus the bound for Ty provides us with 
the estimates needed for the second integral. 

This concludes the proof of Theorem 11.21 □ 


4 Integral estimates 


Proof of Lemma 13.2L We begin by establishing that for G satisfying 

||G||i,l(R) < Cl, llCyllLOOI-R) < C2t ||Cy||il(R) < C^t 

for some positive constants Ci, C 2 , and C3, we have 

p+00 

(4.1) / G{x,t]y)vo{y)dy < GEo{l + 


First, in the event that t is bounded, we have 

^+00 

/ G{x,t-,y)vo{y)dy < ||G||li(] 


< GiEq. 


In the alternative case, for which we take t bounded away from 0, we integrate 
by parts to obtain 

^+00 /*+oo 

/ G{x,t;y)vo{y)dy= Gy{x,t;y)Vo{y)dy, 


where 

Vo{x) := j vo{x)dx, \Vo{y) \ < G^Eq^I + \y\)~^^^. 

J —00 

We have, then, 

1 

Gy{x,t;y)Vo{y)dy 

Gy{x,t-,y)Vo{y)dy 


^+00 


< 


'{\y\<Vi} 


'{\y\>Vi} 


Gy{x,t;y)Vo{y)dy 


< Csf / \Vo{y)\dy + G^Eoil + Vt) ' / \Gy{x,t]y)\dy 

J{\y\<Vi} J{\y\>Vt} 

< G,E,t-^/\ 
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establishing gH). 

It now follows from (EH) that if there exists some ^ 0 so that \x — 
a^t\ < then < C\x — and the hrst estimate is apparent. 

In the event that \x — a^t\ > y/i for all a^, we consider the case in which 
there exists some ^ 0 so that \x — aft\ < et, where e > 0 is sufficiently 
small so that for all j k, \x — a^t\ > rjt, for some rj > 0. In this case, we 
compute 


(4.2) 


'•+CXD 


G{x,t;y)vo{y)dy 


'{\y\< 


\x-aft\ 


G{x,t;y)vo{y)dy + 


} 


'{\y\> 


\x-aftl 


G{x,t;y)vo{y)dy, 


} 


where N will be chosen sufficiently large in the analysis. For the second of 
these last two integrals, integration of G immediately gives an estimate by 

GEoil + \x-aft\)-^/\ 


which is sufficient for our hrst estimate since we are in the case \x—aft\ > y/t. 
Alternatively, for the hrst integral in (EH, we integrate by parts to obtain 
(4.3) 




G{x,t-,y)vo{y)dy 


-} 


< 


x — aft^ x — aft^ 
G{x,t,^^)Vo{- ^ 


N 


N 


+ 


x — aft^ , x — aft^ 
G{x,t, -^)ld)(- 


N 


N 


/ |.-a±t| \Gy{x,t;y)\\Vo{y)\dy. 


It remains to estimate each of these last three terms for each summand in 
the estimates on |G| and \Gy\ from Lemma (2.dl As each case is similar, we 
proceed only with estimates on the convection terms. 

In the case |/ < 0 and for the convection Green’s function estimate, we 
have 


G{x^ t] 


\x - t\ 
N 


)^o(- 


\x - ttj t\ 
N 


< GEo^r^i 


\x — a - t| _ ^ 


(1 + \x 


k=l 
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In the event that j = fc, we have 


I _ I .00 Qju t. . _ ,, 1, 

k - > (x - a^t){l - —), 

which for N sufficiently large leads immediately to an estimate by 

1 /r» 1 /n 

« (1 + |x - , 


for some L sufficiently large. According to the boundedness of 


(4.4) 


X 


~t\ 

' ' Lt 


fV2 


we obtain the claimed estimate. On the other hand, if j ^ k, then |x —a“t| < 
et and \x — a'^t\ > rjt, so that for N sufficiently large, we have exponential 
decay in t. 

We next consider the integral in (D, for which in the case of the con¬ 
vection estimate and for y < 0, we have 



lx —a, . 11 
' J ' 

N 



|x —a . t| 

- rr -“fc *) 


-TJX^ 


Vo{y)\dy. 


Proceeding as in the boundary case, we see that in the case j ^ k, the kernel 
decays at exponential rate in t, while for j = fc, upon integration of Vo{y), 
we have an estimate. 


C'Eor^e-^^^(l + |x - a-ty/2 
< CEo\x-a-t\-^/^. 


For the second estimate in Lemma na we observe from Remark 12.61 the 
estimate 

\dyei{y,t)\ < Ot y] e Mt . 


Integrating by parts, then, we have 

r+oo 

/ ei{y,t)vo{y)dy 


'•-l-CXD 


dyei{y,t)Vo{y)dy 


— OO 


e Mt Eo{l + \y\) ^l^dy 

_ ^ J — CO 
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In this last integrand, we observe the inequality 

(y + alt)'^ , (y + aZ-t)'^ , 

e Mt £'q(1 + |j/|) /^<Ce Mt ^ 0 ( 1 +^) 


from which the claimed estimate is immediate. 

The hnal estimate in Lemma l!01 is proven similarly as the second estimate 
in Lemma 2 of [HZj. □ 

Proof of Lemma 13.3L Looking at ()2.9j) . we notice that in order to estimate 
j +00 jt suffices to estimate 


^ + 00 


n-(x)0{e-’^‘)S,_^.t{-v)Cf(y)vo{y)dy 

< CEoe~'^°^vo{a*t — x) 
<CEQe-^°\l + \a^jt-x\)-i 

<C'Eoe-’'°*(l + |a;|)-t(l + |a;t|)i 

< CEoe~^{l + |x |)“5 


Here we used the crude inequality 

1 


< 


1±M 

1 -|- Icil 


1 + |a + 6| 

and the fact that a* TZ* and are bounded. This gives us (ITT7I) . Estimate 
(ITTHl) is obtained similarly. □ 

Proof of Lemma I3.4L For Lemma |S31 the proof of each estimate requires 
the analysis of several cases. We proceed by carrying out detailed calculations 
in the most delicate cases and sufficing to indicate the appropriate arguments 
in the others. In particular, we will always consider the case x,y < 0. The 
case 1 / < 0 < a: is similar (though certainly not identical) to the reflection 
case for x,y < 0, and the estimates for y > 0 are entirely symmetric. The 
analysis of each type of kernel—convection, transmission, and scattering— 
is similar, and we carry out details only in the case of convection. (This 
terminology is reviewed in Remark |2.4j) . 

Nonlinearity (1 + We begin by estimating integrals, 

(x-y-a~ (t-s))^ 


(4.6) 



(t-s) 


-V 


0 J — |s 


M{t-s) s|) 3/4g i/2('x_|_5) ^/^dyds. 
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In the event that |x| > \ai |t, we write 


(4.7) X — y — ttj (t — s) = {x — ttit) — {y — ai s) + (a;^ — )(t — s), 

and observe that in the current setting {x < 0, y E [—|a7|s,0], < a^), 

there is no cancellation between these three summands. Integrating (1 + ||/ — 
for s G [0,t/2] and integrating the kernel for s G [t/2,f], we obtain 
an estimate by 

Jo 

(4-8) + f {t - s)-^/^ds 

Jtl2 

(x-a~ t )2 

< Ct ^ ln(e + t)e ^ 


which is sufficient by the argument of (jOl). We observe that the seeming 
blowup as f —0 can be eliminated. Integrating the kernel in (gH), we have 
an estimate by 

C f (t — s)~^^'^s~^^‘^ds 

fA n\ Jo 

s~^^'^ds + C2t~^^‘^ f ~ s)~^^‘^ds < C. 

Jtl2 

For |x| < |a7|t, we write 

(4.10) X - y - aj{t- s) = {x - a~{t - s) - a^s) - {y - a^s), 



from which we observe the inequality 


{x-y-a^ {t-B)y 


(4.11) 


e M(t-a) s|) 3/4 

< c 


{x-a. b) 


(x-y-aj (t-B))‘^ 

M(t-B) (^x + I// — a^s|)“^/^e~'^ 


M(t-B) 


(x-y-a. (t-B))^ 


+ e (1 + ly - a;;.s| + |x - a-(t - s) - a;;.s|) 


for some M sufficiently large and e > 0 sufficiently small. The analysis can 
now be divided into three cases: (i) < 0 < aj, (ii) ^ ^^*4 
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(iii) ttj < CLf. <0. Each critical argument appears in Case (ii), and so we 
consider only it in detail. For the hrst estimate in dnn), we have 


t rO 



(x-y-a. (t-s))^ (a:-a . (t-s)-a, a)'^ 


{t — s) M(t-S) g M(t-s) 

0 ^ —|aj”|s 

X (1 + ly — + s)~^^‘^dyds. 


In the event that |x| > 1% |t, we write 


(4.12) X - ttj (t - s) - a^s = {x - a^t) - {ttj -a^)(t-s), 

for which there is no cancellation between the summands, and we can proceed 
exactly as in (USD- On the other hand, for |a:| < \t, we write 

(4.13) X — aj{t — s) — a^s = {x — ajt) — — a~)s, 


for which again there is no cancellation between summands, and we can 
proceed exactly as in (gSl). For the critical case |a^- \t < |a;| < \t, we 

subdivide the analysis into cases: s G [0, t/2] and s G [t/2, t]. For s G [0, t/2], 
we observe through iH the estimate 


(4.14) 


(x-aj (t-s)-aj^ s)^ 

M(t-s) 


(1 + ^) 


-3/4 


< c 


(x-aj 


(x-a. {t-s)-a^ s)^ 


e M(t-s) (1 + s) e Mit-s) + |a: — a ■ t|) 


-d 1-3/4 


For the hrst, we can proceed exactly as in (gSl), while for the second, we 
integrate (1 + |?/ — % sj)"^/'^ to obtain the estimate 


Cit (1 + |x - aj t\) 


-dl-3/4 


rt/2 


(x-a. (t-s)-a^ s)^ 


+ sf/^ds 


< Ct ^'^^(l + |x — ajt|) 


where in establishing this last inequality we have integrated 

(x-a~ (t-s)-a~ s)2 


in s. For s G [t/2,t], we observe through ()4.12|1 the estimate 
(4.15) 

(x-a. (t-s)-a, s)2 

{t - 


< C 


(x-a, t) 


(x — a^ (t —s) —s)^ 


(t — S) 3/4g Mit-s) 
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For the first of these estimates, we proceed exactly as in dUD, while for the 
second, we integrate the kernel in y to obtain an estimate by 




C2\x-a,^t\ M(t-s) ds 

it/2 


The apparent blowup as x approaches t can be removed similarly as in 

(ED. 

For the second estimate in (gUD, we have 


t rO 



(x-y-a. (t-s))^ 


{t — S) M(t-a) 

^0 ^ —|a]”|s 

X (1 + I// - a^s| + |x - a“(t - s) 


«fcS|) ^/^(1 + s) ^/‘^dyds. 


Again, we proceed in detail, only in Case (ii), a'^ < a~ < 0. In the event 
that |a;| > |a^|t (though keeping in mind that we have already considered the 
case |x| > |ai“|t), we use (I4.12|l . for which there is no cancellation between 
the summands, and upon integration of the kernel, we have an estimate by 

C{l + \x - [ (t - + s)“^/^cis 

Jo 

.tl2 

(4.16) <C^{\^\x-alt\)-^lH-^l‘^ \ s-i/2(l + s)-3/4ds 

Jo 

+6*2(1 + |a; - [ {t - s)~^^‘^ds, 

Jt/2 


from which we obtain an estimate by 'ipi. Likewise, in the case |a:| < \aj \t, we 
use (inm . for which there is no cancellation between summands, and we can 
proceed similarly as in (jmni). For the critical case |a^- \t < |a:| < |a^ |t, we 
subdivide the analysis into cases: s G [0,f/2] and s G [t/2,t]. For s G [0,f/2], 
we observe through the estimate 


(1 + I// - s| + |a; - a • {t - s) - s|) ^^'^(1 + s) 


< C 


(1 + \x-aA\) ^/^(1 + s) 


+ {1 + \y - s\ + \x - ttj {t - s) - s|) ^^^(1 + s + \x - t\) 
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For the first of these estimates, we can proceed as in (iron . while for the 
second, we integrate the kernel in y to obtain an estimate by 

/•*/2 

Ci{l + \x — / (1 + |a: — a~{t — s) — 

Jo 

< C{1 + + lx — 

For s G [t/2,t], we observe throngh (I4.12|l the estimate 
{t - + \y - a^s\ + \x - aj{t - s) - 

+ \x - + \y - a;^s| + \x - aJ{t - s) - 

For the hrst of these estimates, we proceed exactly as in (UHl), while for the 
second, we integrate the kernel in y to obtain an estimate by 

C 2 \x - [ {t - + |x - a“(t - s) - 

Jt/2 

< C\x - 

Nonlinearity (1 + We next consider integrals of the form 

(4.17)^ 

(t — s)~^e M{t-s) (X _|_ _ Q-_ 5 | -|_ s^/2^-3/2^-i/2j']^ gY^^'^dyds. 



0 j —OQ 


In this case, we again observe (ginD, from which we obtain the estimate 
(4.18) 

{x-y-a. (t-s))^ 

e M(t-s) (^x + li/ — a^sj + sh2^-3/2 


< c 


(x-a. (t-s)-a, sy (a,_y_a 


M(t-s) g ^ M(t-s) (^x + |j/— s| + 


{x-y-a- {t-a)Y 


+ e (1 + li/- a^s| + |x - a -(t - s) - a^s| + s / ) 


l/2\-3/2 


As in onr analysis of the nonlinearity (1 + s) we focns on the case 

a'^ < a~ < 0, compared to which the other cases are either similar or more 
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straightforward. For the hrst estimate in (ICTl . we have 


t rO ^ (x-y-a. (t-s))2 (x-a^ {t-B)-a^ a)^ 

(t— M(t-s) g M(t-B) 

'0 J —oo 

X (1 + |j/ — a^s| + + s)~^^'^dyds. 



In the event that |x| > |a^ |t, we employ (14.1211 . for which there is no cancel¬ 
lation between the summands. In this case, we have an estimate by 

(x-a-t)2 /-Fa , , , 

Jo 

(4.19) + ^*2(1 + 1/2(1 + r {t - s)-i/2ds 

Jt/2 

{x — aTt')^ 

< Ct~^ ln(e -1- t)e « , 


where the seeming blowup as t —> 0 can be eliminated as in (USD. In the 
case |x| < |a~|t, we employ (j4.1 ,‘fjl for which again there is no cancellation 
between the summands and we can proceed as in (ICTll . 

We turn next to the critical case |aj|t < |x| < |a^|t, for which we divide 
the analysis into cases s G [0, t/2] and s G [t/2, t]. For s G [0, t/2], we observe 
through ()4.13|1 the inequality 


{x-a. (t-s)-aj. s)2 


-s-1/2(1 + 5)-1/4 


< c 

+ e 


(x-a. t) 




Lt S 


- 1/2 


(1+^) 


-1/4 


{x-a. (t-s)-a^ s)2 


|x —a-t| 1/^(1-I-|x — a-1|) 1/^ 


For the hrst of these last two estimates, we can proceed as in (j4.19|l . while 
for the second, we have an estimate of the form 


Cit i|x —a-f| 1/^(1-|-|a: — a• t|) 


-hn-1/4 


l•t/2 (x-a. (t-s)-a, s)'^ 


-(l + si/2)-i/2ds 


< Ct i/^|x —ajf| 1/^(1-I-|a: — ajt|) i/^^, 
where as usual the singular behavior can be removed. For s G [t/2,t], we 
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observe through the inequality 




(t-s) 




M(t-B) 


< c 


(x-a. (t-s)-aj^ s)2 


\X — aj^t\ M{t-s) (^f — g'j 3/4g 


(a;-a^ t)^ 


For the second of these last two estimates, we can proceed as in (ICTll . while 
for the hrst we have an estimate of the form 


{x-a. s)2 


Cilx - + ds 

Jt/2 

< C\x - 


For the second estimate in (jd.lHj) . we have 


t rO ^ (x-y-a. (t-a)y 

{t — S)~^e M(t-B) 


Jo J —oo 

X (1 + |j/ — a^s| + |x — aj{t — s) — a^s| + + s)~^^'^dyds. 

Focusing on the case < a~ < 0, we observe as before that for |a:| > |a^|t, 
we have no cancellation between the summands in ()4.12|1 . and consequently 
there holds 


(1 + I?/ - s| + |x - Oj {t - s) - s\ + 

< (7(1 + \y - a^s\ + \x - a'^t\ + (t - s) + 

from which, upon integration of the kernel, we have an estimate by 

C.(l + |x - alt\ + l\t - s)-'/2s-i/2(l + s)-'/* 

< (7(1 + t)-^/^(l + |x - a^t\ + t^/2)-3/2_ 

In the event that |a;| < |a~|t, we have no cancellation between the summands 
in (mi, and consequently there holds 

(1 + ly - a^s| + \x - aj{t- s) - a^s| + sh2^-3/2 
< (7(1 + \y - a^sl + |a; - 
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from which we have an estimate by 


Cit ^(1 + \x - a - 1\) 


-fh-V2 


ft/2 




(4.21) 


+ C 2 t ^/^(1 + t) + |x - a- 1| + t) / (t-s) 

Jtl2 

<Ct-^/\l + \x-a-t\)-^/\ 


In the critical case \t < |x| < |a^ |t, we snbdivide the analysis farther into 
cases s G [0,t/2] and s G [t/2,t]. In the case s G [0,t/2], we observe throngh 
the estimate 

(1 + ll/ - a^s| + |x - a~{t - s) - a'^s\ + 


< C 


(1 + |r/ — s| + |x — a • {t — s) — s| + |x — a ■ t\ + 


+ (1 + |r/ — Ofc s| + |x — a“(t — s) — a^sl + 5^2^-372 


X (s + |x — a ■ t\) + s + |x — a• t|) 


- A-th-V4 


For the hrst of these last two estimates, we can proceed as in don, while 
for the second, we obtain an estimate by 


Cit ^|x —a-tl + |x — a■ t|) 


j 

ft/2 


X 


(1 + |x — ttj {t — s) — ttf, s\ + 

Jo 

< Ct~^^‘^\x — + |x — 

where the apparent singnlarities can be eliminated as in (jOl). In the case 
s G [t/2,t], we observe through (j4.12ll the estimate 

(t - + \y - a^s| + |x - a“(t - s) - a^s| + 5^2^-372 


< C 




ttk s 


+ |x - a ■ (t-s) - ttf, s\ + s^/^) 


+ (t - s) + \y - ttk s\ + \x - a^t + s^/^) 


For the second of the last two estimates, keeping in mind that here s G [t/2, t], 
we can proceed as in (ICTD . while for the hrst we integrate the kernel to 
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obtain an estimate by 

Ci\x - Ofc 

X [ {t — + \x — a'^{t — s) — a^s| + |x — a'^t\ + 

Jt/2 

< c\x - 

where the apparent singnlarities can be eliminated as in (USD- 

Nonlinearity (1 + We next consider integrals of the form 

(4.22) r r (f - + \y\)-^/^dyds. 

Jo J —\a-^^ \s 

In the case |a:| > |ah|t, we observe that there is no cancellation between 
snmmands in decomposition (EH), and conseqnently, we have an estimate by 


(4.23) 


(x —a, t)^ 


1 1 /r, (^-a, t) 

+ 6 * 2(1 + ~ 


2 /•t 


hl2 


{t — s) ^^‘^ds 




< Ct ^ ln(e + t)e ^ 


where the apparent blownp as f —> 0 can be removed as in (ED. For the 
remainder of the analysis, we restrict onr attention to the case |a;| < |ahl^ 
and aj < 0 (the case a“ > 0 is similar and more direct). Observing that for 
y G [—|ah|s,0], we have 

(1 + s) ^ < ^(l + s + li/l) \ 


we have the ineqnality 


(x-y-a. (t-s))2 

e M{t-s) + s)“^s“^^^(l + \y\)~^^^ 


{j/e[-|aj p,0]} 


(4.24) 


< C 


(x-y-a^ (t-s)) (*-»)) 


M(t-s) g Mit-s) (1 + S) ^/^(l + |j/|) 


(x-y-a. (t-a))^ 


+ e (1 + S + ll/l + |X - Oj (f - s)|) 

X (s+ \y\ + \x - aUt - s)\)-^/^{l + \y\ + \x - aJ{t - s)\)-^/^ 
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For the first estimate in we observe that for |a:| > \aj \t, there is no 

cancellation between x — a~t and a“s, and consequently we have an estimate 
by 

C,t-^e -^ / {l + s)-^/^s-^/^ds 

Jo 

(4-25) +c,(l + r (t - 

Jt/2 

< Ct~^ ln(e + t)e , 


where the apparent blowup as t ^ 0 can be removed as in (jOj). In the case 
|a;| < |a“|t, we divide the analysis into the cases, s G [0,t/2] and s G [t/2,t]. 
For s G [0,t/2], we have the inequality 


(x-a^ (t-s)) 


M{t- 


— (l + s)-V^/2 


< C 


(x-a^ t) 


(x:-a^ (t-s)) 


Lt e M(t-s) + s) 


(x-aj (t-s)) 


+ e M(t-s) -|_ 5 I 2 ; _ Q + 1 ) jf. — a-1\) 




For the hrst of these last two estimates, we proceed as in while for 

the second, upon integration of (1 + \y\)~^y, we have an estimate on ()4.22j) 
by 

ptj‘2 {x — a . (£ —s))^ 

+ |x — / e M(t-s) 

Jo 

X (1 + s + |x — + |x — + s^^'^ds 

< + |x - a“£|)“^/^. 

In the case s G [£/2,£], we observe the inequality 

(^ _ 5)-l/2g M(t-s) <C\x\~^^‘^e M(t-s) ^ 


from which, upon integration of the kernel, we obtain an estimate on 

by 


C'2(l+ £)-'£ 




{x — aj{t — s))^ 


M(t-s) dg 


't/2 


<Ct-^\x\-^/\ 
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For the second estimate in (j4.24|l . we observe that in the case |x| > \aj\t, 
we have no cancellation between {x — ajt) and ajs, and consequently can 
estimate 


(4.26) 





+C 2 t + \x - a,t\) [ {t - s) 

Jt/2 


< C{1 + |x - aj t\) 


1/2 


which suffices since the case |x —a“t| < C^/i requires only decay. In the 
case |a;| < \aj\t, we subdivide the analysis further into the cases s G [0,t/2] 
and s G [t/2,t]. For s G [0,t/2], we observe the inequality 


(1 + s + |x - Qj (t - s)\) ^ < (7(1 + s + |x - a^-1|) \ 


from which, upon integration of the kernel, we obtain an estimate by 


rt/2 


Cit ^^^(l+|x—a-1|) ^ / s ^/^(l+|x—a■ (t—s)|)(is < (7f ^^^(l+|x—a■ t|) ^ 


For s G [t/2,t], we observe the inequality 


{t - s) ^/^(l + \y\ + |x - a- (f - s)|) 




< 


x\ ^/^(l + |i/| + |x - a-(t - s)|) + (t - s) ^/^(l + |i/| + |x|) 


For the hrst of these last two estimates, we obtain an estimate on (j4.22j] by 
(72(1 + t)-H-^/^\x\-^/^ f (1 + |x - a-{t - s)\)-^/‘^ds < C{1 + t)-^\x\-^/\ 

Jt/2 

while for the second, we have an estimate by 

^2(1 + + |a;|)-V2 < ^(l + t)-^\x\-^/^. 

Jt/2 


Nonlinearity s ^/^(l + s) We next consider integrals of the form 


(4.27) 



{t-s) 




{x-y-a. {t-s)y 


M{t-s) 


■(1 + ^) 


-■i/A-l/2- 




dyds. 


0 J — oo 
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In this case we observe from Lemma 6 of [HZ 


the equality 


(x-y-a. (t-s))2 (y-a- b)^ 

g M(t-s) g Ms 

(4.28) 

(x-a. (t-3)-a^sy £ -^2 

— g Mt g A'/s(t —s) 


from which direct integration over y leads to an estimate by 

(x — a ■ (t — s) — a, sY 

Ct-Y^ / - s)-Y\l + s)-Y^e - ^ M. 

Jo 


As in the previous analyses, we focus on the case < aj < 0, and note 
that analysis of the remaining case a“ < < 0 is similar. In the event that 

> lofc |£, we observe that there is no cancellation between the summands 
of ()4.12j) and we have an estimate by 


Cit Yt 


(x-a, tf H/S 


(1 + s) Y^ds 


(4.29) 


(x-a-). Y 


+ C 2 t-Y\l + t)-Y\ -/ {t-s)-Y^ds 

Jt/2 


< Ct-Y^e - 


For |a;| < |a^|£, we divide the analysis into cases, s G [0,£/2] and s G [£/2,£]. 
For s G [0,£/2], use (I4.18|l . for which the summands cancel and we have the 
estimate (gm. For the hrst estimate in (gm, we can proceed exactly as 
in while for the second, we have an estimate by 


Cit ^(1 + |a; - a - £|) 


ft/2 


(x-a. (t-a)-a s)^ 


' ds 


< Ct -\-\x — ajt\) 


where the apparent blowup as £ —>■ 0 can be removed as in gSl). For s G 
[£/2,£], we observe through (I4.12|l the estimate ()4.15j) . For the hrst estimate 
in glSD, we proceed as in g^ni), while for the second we have an estimate 
by 





(x-a. (t-B)-a, b)^ 

-- ds 
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Nonlinearity (1 + s) We next consider integrals of the form 

pt nO (^x—y — a ■ (t — s))^ 

(4.30) / / {t — s)~^e M(t-s) + s)~^dyds. 

Jo J —oo 

First, we observe that in the event \y\ > |a]”|s, we have exponential decay 
in both y and s, from which the claimed estimate readily follows. For what 
remains, then, we focus on the case \y\ < |a]”|s. In the event that |a;| > |a]”|t, 
we write 


X - y - ttj {t - s) = {x - a^t) - {y - Oi s) - (a^ - a^){t - s), 

for which there is no cancellation between the summands (in the case y G 
[a]”s,0]), and we immediately arrive at an estimate by 


In the case |a:| < |t, we focus on the critical case < 0. Here, we observe 

the estimate 


{x-y-a (t-sN 


(4.31) 


M(t-a) 


^-vlvl 


< c 


(x-a^ (t-s)) 
M(t-s) 




for some constants M > 0, rji > 0, and '>]2 > 0. For the hrst estimate in 
(lOTl) . we have 




(t-s) 




(x-a^ (t-s))2 

+ s)~^dyds. 


In the event that |a;| > |a^. \t, we write 

(4.32) X - aj{t - s) = {x - ajt) + a“s, 

for which there is no cancellation between the summands, and we arrive at 
an estimate by 

(4.33) 

(x-a~t)^ ri/2 (x-a~t)'^ 

Cit~^e / {1 + s)~^ds + € 2(1 + t)~^e / {t — s)~^ds 

Jo Jt/2 

+ C 3 (l+t)“^e / {t — s)~^^‘^ds < Ct~^\n.{e + t)e . 

Jt-i 
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For |a:| < |a^- \t, we divide the analysis into cases, s G [0,t/2] and s G [t/2,t]. 
In the case s G [0,t/2], we observe the cancellation between snmmands in 
(lO^ . which leads to the estimate 


(x-a^ {t-s)Y 


e M(t-s) + s) 


-1 


< C 


(x-a^ t) 


(x-a^ (t-s)y 


Lt (1 + S) ^ + e M(t-s) (^x + \x 



For the hrst of these last two estimates, we proceed exactly as in flOHll . while 
for the second, we have an estimate by 

/■i/2 (x-a~(t-s)Y 

+ |x — ajt\)~^ / e M(t-s) f/g < + |x — 

Jo 

For s G [t/2,t], we observe that for |a“|(f — s) < (l/2)|a:|, we have the 
estimate 

Y-a~(t-s)Y ^ ^2 (2:-aT(t-s))2 

(/ _ g)-l/2g M(t-s) < _ 5)-l/2g-rRg M(t-s) ^ 


while for \aj |(f — s) > (l/2)|x|, we have the estimate 


{x-a^ (t-s)y 


(x-a^ (t-s)Y 


{t — s) M(t-a) < C\x\ M(t-s) 

For the second of these last two estimates, we obtain an estimate 


{x-a^ (t-B)Y 


it/2 

For the hrst, we have the same decay in t, with exponential decay in |x|. For 
the second estimate in ()4.dl|l . we have 


t /-O 


(t - + s)-^dyds. 


Jo J — oo 

In the event that |a:| > \a~\t, we observe that there is no cancellation between 
the snmmands of (0321), so that 

^-r^l\x-a~ {t-s)\ ^ (J^-y4,\a~\s 
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which leads immediately to an estimate better than 'ip 2 - For |a:| < |a^- \t, we 
proceed similarly as with the hrst estimate in (lOTD . 

This concludes our analysis of the main case, the hrst estimate in Lemma 

in 

Second estimate of Lemma \S.4\ For the second estimate in Lemma Id.41 
we consider in detail the case of nonlinearity dy{(p^)‘^, for which the other 
cases are similar. Following the analysis of |H,aj . we divide the integration 
over s as 


(4.34) 

pt /•+00 ^ py/t ^+00 

/ / G{x,t-s;y)dy{ip'^fdyds= / / 

Jo J—oo Jo J—oo 

pt—y/t /•+00 

+ / / G{x,t - s;y)dy{(p^fdyds 

J y/t J —OO 

pt /•+00 

^k\2 


G{x, t — s] y)dy{ip^)‘^dyds 


+ 


’t—y/t J — C 


G{x,t — s-,y)dy{(p ydyds. 


For the hrst integral in (im . we focus on the case y < and on the the 
convection term of the Green’s function, for which we must estimate integrals 
of the form 


o^/t rO 


(t-s) 


(x-y-a. 

-'e-^ 


M(t-s) 


(1 + .) 


-'e- 




dyds. 


According to Lemma 6 from [HZ], we can write (imi . from which direct 
integration over y leads to an estimate by 


(4.36) 


/‘V^ (x — a- (t —s) —a, s)^ 

Gt-^/^ / {t - - ^ 


ds. 


We have three cases to consider: (1) < 0 < a^- , (2) < 0, and 

(3) aj < < 0; here we point out that the case = aj does not arise: 

the reason is that, according to (USD, the term in $ occurs only in 

r~ direction for j ^ k. Hence if we consider the detailed description of S in 
Green function bounds dnsD, we can see that the case does not 

arise (since Ijr^ = 0). We will focus on the case af < a~ < 0, for which the 
remaining cases are similar. First, in the event that \x — aft\ < Gy/t we can 
conclude decay of the form 

(x-a~ t)2 

e 
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by boundedness. In the case \x — t\ > Cy/t, for C sufficiently large, and 
for s G [0, \/i], we have 

\x - a~{t - s) - a'^s\ = |(a: - ajt) - (a^ - aj)s\ > Ci|x - a~t\, 
for which we can estimate (1031) for f > 1 by 

< «-■(! + 

For the third integral in ()4.d4j) . we can proceed similarly to obtain an estimate 
by 

, , {x — a, £)^ 

Ct-^/\l + t)-^/^e - 

In the case of the second integral in ()4.d4|l . we closely follow the approach 
of |Liulj . using the framework of jR,aj . Dehning for some fixed k 4>{x,t) as a 
Burgers kernel, we have 

= -7fc for t > -1, 

0(x + a^, 1) = t =- 1 , 

for which we have (p^(x, t) = 0(x — a'^t, t), where the (p^(x, t) are as in (jl.lbjl . 
In this notation, the second integral in (lOl) becomes 

pt—y/t p+OO 

(4.37) / / G{x,t - s;y){(j){y - a^s,sY)ydyds. 

J \/t J —OO 

The idea of Liu is to take advantage of the observation that time derivatives 
of the heat kernel decay at the same algebraic rate as two space derivatives, 
and also of the relation between time and space derivatives of G and 0. 
Intuitively, we can think that we would like to integrate (1071) by parts in s 
to put time derivatives on G, but in order to facilitate this, we would like the 
derivatives on 0^ to be with respect to s. Though convecting kernels such 
as G do not enjoy this property of fast-decaying time derivatives, we observe 
that under a change of integration variable, the leading order term in G can 
be converted into a heat kernel. More precisely, according to Proposition |2i21 
the leading order contribution to the scattering piece of the Green’s function 
is given hj g{x — aJt, t), where g{x, t) is the heat kernel 

g{x,t;y) -.= , 


(4.36) Git 




(x-a. t)^ 
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for a constant c = {Ij )*“'/y 47r/?j . Under the change of variables ^ = 

y + aj{t — s), the second integral in (I4.d4j) becomes, for this leading order 
term, 


/• + O0 




g{x,t- s;O(0(^ - aj {t - s) - a^s, sf)^dids. 


Following jH,aj . we write for g = g{x,T,^) and 0 = 0(^, r). 


(^g{x, t - s; 0 - S' ^ 

+ («7 “ ^ 0(7>(^ - aj{t-s)- s, s)2)^ 

+ t - s] O(0(^ - aj {t-s)- s, s)^)^, 

from which we have 

(«7 - ak) 9 {x,t - s;^){(l){^ - aj{t - s) - a^s)^)^ 

= (^9{x, t - s] 0 (7>(^ - “7 {t-s)- s, s)^)) ^ 


(4.38) 


+ gr{x,t- s]C){(t){i-aj {t-s) -afcS,s)^) 
- 9{x, t - s; 0(7>(^ -aj{t-s)- s, s)^)^. 


We proceed now by analyzing (nnder integration) each term on the right 
hand side of ()4.38j) . For the hrst, we have 


rt—\/t ^+oo 

' \/I J —oc 


(4.39) 


^H-CXD 


’ —oo 

P + OO 


(^g{x,t- s;^){(p{^- ttj {t-s) - ai^s,s)^)^ d^ds 
g{x, Vt; O(0(^ -ajVi-a^{t-Vi),t- Vi)‘^)d^ 


+ / 9{x, t - Vi] O(0(^ - a ■ {t-Vi)- Vi, Vi)'^)dC 


For the hrst integral in (I4.39j) . we have an estimate by 

(i-a~ Vi-a7 (t-Vi))'^ 


f + oo 


C I {Vi) (1 + (t — Vi)) M(t-V) 
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Writing now ^ = z + \/t, we can rewrite this last integral as 


C 



{x — z — a ■ >/£)^ 

"'/i (1 + (t — \/t))~^e 


M(t-Vi) 


We have, then, according to Lemma 6 of [HZj . with t — s replaced by \/t, 

(x-z-aj (z-a^(t-\/t))^ 

g Mv^ g 

(x-a~Vt-a~(t-y/t))'^ _ f +11^, )x/t(t-Vr)^ 

= g Wt g MVt(t—Vt) t ^ 


Upon integration in then, we have an estimate by 

{x — a~ \/t — a 7 {t — y/typ' 

C<rV2(i + (t _ - x/t)V2e--. 

In the event that |x — a^t| < CaA, we have exponential decay 

(a:-art)^ 

e Mt 

from boundedness, while for |x — a^t| > C\/i, C sufficiently large, we have 
\x - ajVi - a^{t - \/i)\ = |(x - a^t) - {aj - a^)\/i\ > c|x - a^t\. 

In either case, we conclude for t sufficiently large an estimate 


The second integral in (jd.dbj) can be analyzed similarly, and we obtain an 
estimate of the same form. For the second integral arising from , we 

have 


p+oo 


fvt 


{t — s) + s) ^e' 


(5-0,- (t-s)-Oi. s) 


'd^ds. 


Returning to our original coordinates ^ = y + a, {t — s), we have 


/‘+00 


(x-y-a. (t-s))^ 


'Vt 






d^ds. 
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Proceeding now exactly as in the derivation of , we arrive at an estimate 
by 


(4.40) 






(x-a. (t-s)-a, sy 
^-U1 , - ^ Mt 


'Vi 




ds. 


As in previons argnments, we will focns on the case < a“ < 0, for 
which the remaining cases are similar. First, in the event that |x| > \a^\t, 
we observe that there is no cancellation between snmmands in (gH, and 
conseqnently that we can estimate (lonii for t large enongh so that \/t <t/2 
by 


(x-a, t)2 


-/ {l + s)-^/^ds 


'Vi 


(4.41) 


r^-Vi 


{t — s) 


lt/2 


(x-a. (t-s)-a^ s)2 


'ds 


< cr^e 




For |a;| < |a~|t, we observe that there is no cancellation between snmmands 
in (I4.1d|) . and conseqnently we can estimate (I4.4()|l similarly as in (I4.41|l to 
obtain an estimate by 


For the critical case \aj \t < |x| < |a^ \t, we divide the analysis into the cases 
s G [Vt,t/2] and s G [t/2,t — Vi]. For s G [\/i,t/2], we observe throngh 
(imni the ineqnality 


(x-a . (t-s)-a^ s)2 


< c 


{x-a. 4)2 


(1+ + e- 


(x-a . (t-s)-a^ s)2 


(1 + \x 



For the hrst of these last two estimates, we proceed as in (lOTl) . while for 
the second, we have an estimate by 

1*^/2 {x — a~{t — s)—aTs)^ 

t|)-V2 / e- ^ m4 ^ ds<Ct-^/\l+tY/^ 

JVi 


(l+lx-a^. t|) 


Cit (l+lx-Oj 


- 1/2 
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For s G [t/2, t — \/i], we observe through (iroi) the estimate 




(t-s) 


(x-a. (t-s)-a^ s)2 


< c 


\x — a^t\ 


(x-a. (t-s)-aj, s)2 


+ (t-s)-'/V 


(x-a, 4)2 


For the second of these last two estimates we can proceed as in (14.4111 to 
determine an estimate by 

Ct~^e eF 

while for the hrst, we have an estimate by 


(x-a^ t)2 






(x-a^ (t-a)-a^ s)2 


ds 


'Vt 


< C{1 + t)-^/^\x - 


For the hnal expression on the right hand side of ()4.d8jl . we have integrals of 
the form 

pt—\/t r+OQ 


Wt 


— (^ —a - (t —s) —a, s)'^ 

{t — (1 + s)~‘^e ' Ms d^ds. 


Re-writing in our original variable ^ = y + aAt — s)^we have 


f+OO 


(x-y-a. (t-a)Y 


'Vi 


[t — s) Mit-s) -1- s) 




d^ds. 


Proceeding now exactly as in the derivation of (j4.4()jl . we arrive at an estimate 
by 


(4.42) 


pt y/t (x — a- (£ —s) —a, s)^ 

/ (1 + - ^ M4 ^ 

Jvi 


ds. 


For the cases |x| > |a^|£ and |a;| < |aj|t, we can proceed as in (jd.dljl to 
arrive at estimates of the form while in the case |a“|f < |a;| < |a^|£, 

we divide the analysis into case, s G [\/£,£/2] and s G [£/2,£ — -y/t]. For 
SG [v^,f/2], we observe through (imH) the estimate 


(x-a^ (t-a)-a^ s)2 


(l+«) 


-3/2 


< c 


(x-a. t)2 


it (1 + s) + {l + \x - aA\+t^^‘^) 


(x-°-j (t-s)-aj. s)2 
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For the first of these last two estimates, we proceed as in dOTl) . while for 
the second we have an estimate by 


/•h2 

JVi 

< C{1 + |x - ajtl 


{x-a. (t-s)-a, s)^ 

-- ds 


For the case s G [t/2,t — \/t], we similarly arrive at an estimate by (1 + 

Excited estimates. We turn now to the estimates in Lemma (D that 
involve the excited terms ei{y,t). Observing that in the Lax and overcom- 
pressive cases, we have the estimate 

\dyei{y,t)\ < Ot ^ e Mt , 
a->o 


integrals in the third estimate of Lemma Id.41 take the form 




(t 


, f.y+°'k 

_s)-V2g - M(t-S) 


-bi |s 


^{s,y)dyds, 


which differ from those of the previous analysis only by factor (t — s) 
Proceeding almost exactly as in our analysis of the hrst estimate of Lemma 
HI we determine an estimate on these integrals of 0(1 + t) 

Excited diffusion wave estimates. The third and fourth estimates in 
Lemma El are similar and straightforward to prove, and we consider only 
the third. As each summand in <F(|/, s) can be dealt with similarly (see (Id.711 - 
(ESD), we focus on the nonlinearity for which integration by parts 

in y and the estimates of Remark 12.hi yield integrals of the form 



(y + a. 

1/2^-^ 


M{t-s 


(1 + ^) 


-'e- 


(H-a, s)^ 


'dyds, 


with aj > 0 (the excited terms correspond with mass convecting into the 
shock layer) and a?; < 0 (diffusion waves care mass away from the shock 
layer). According to Lemma 6 of [HZ], we have the relation 

a — ^Y) (y —(t — + s)^ 

g M(t — s) g Ms = g Wt 


g Ms{t — s) 


f , +0-1. )(t-s)s. o 

- 1 - ? 
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from which integration over y leads to an estimate of the form 




ds. 


In either the case s G [0,^/7] or s G [t — t/'y,t], for 7 sufficiently large, we 
have 

\a~{t — s) + a'^s\ > r]t, 

for some 7 > 0, through which we have exponential decay in t. In the case 
s G [t/'y,t — t/'y], we integrate the kernel in s to obtain an estimate by 
C(l + t)-V2. 

For the hnal estimate in Lemma na we have for |/ < 0 integrals of the 
form 

rt rO (y+a-(t-s)f 

/ / (t — s)~^e ^(i/^s)dyds, 

Jo J —00 

where < 0. In the case of nonlinearities (1 + + '02 + 

we can proceed by setting a; = 0 in the scattering estimates to obtain an 
estimate by + In the case of nonlinearity (1 + + \y\)~^^‘^, 

we have the estimate 

(v + ajy (t-s))^ (v + a~ (t-s))^ 

e M(t-s) _|_ |y|)-i/2 < ( 7 g M(t-s) (^x + |r/| + (t — 
from which we immediately obtain an estimate by 

Cit-1/2(X + t)-V2 + s)-^S-^/^ds 


+ C'2(l + {t- + it- s))-^/^ds 

Jt/2 

< C{l + t)-\ 

Finally, in the case of nonlinarity (1 + we have the estimate 

{y-\-aT (£ —s))^ 

g- M(t-s) — Q-v\y\ < (7g-’?i(i-«)g-»?2hl^ 

from which we immediately obtain an estimate by 

pt/2 ft 

(Fie-^* / e-^(*-")ds + C'2(l + £)“W 

Jo Jt/2 


<C{l+t) 


-1 
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This concludes the proof of Lemma Id .41 □ 

Proof of Lemma 13 . 5 L To show (jd .21111 we need to estimate 

pt r + CX) 

Jo J-cx 

<C [ (e-^“(*-"))|/(-a: + a;(t-s),s)|ds 
Jo 

by the boundedness of 71 * and £**. A typical term of 'ipi is a term of the 
form (1 + + |x — a^f|)“2 (or with replaced by af). Hence for the 

terms coming from 'ipi the above is of the order 

(X -j_ s)“ 2 (i + |x — a*{t — s) — a~s\)~^ds 

w(i-s) 5 -|(x s^-|(x 1 ^; _ _ (a* _ al){t - s)\)~^ds-, 

now, using (EH), this is smaller than 

g- 7 ?od-s)g-^(-X ^ s)“^(l + |x - a“t|)“5(|l + (a* - a7){t - s)|)5cis. 

Notice that (a* — a~) < C . Now we use half of rj to neutralize t — s, and to 
get 

pt 

(1 + |a; — + s)~"^ds 

Jo 

< (1 + |a; - a“t|)“t(i + t)-\ 

obviously absorbable in -01 +■02 + a. For the terms coming from '02, we notice 
that, by inequality a? + h’^ > 2ab, 

(1 + ly - a~s\ + \/s)“t < C{1 + s)“^(l + \y - a 0 s|)“L 

Hence applying the same calculations as for ipi gives us an estimate of (1 + 
\x — a0t|)“4(l + t)“8. For a, we apply the same procedure, and using again 
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(03, we have an estimate by 


^ I ^ 2(1 + s) 4(1 + |a; — a*(t — s)|) 2 ds 

pt 

<C / + s)“3(1 + |x|)“^(l + |a*(t — s)|)^(is 

Jo 

<C(l + |a ;|)“5 f + s)~^ds 

Jo 

<c{i + \x\)-'Hi + t)-\ 

This hnishes the proof of ()3.26j) and ()3.27|) . The proof of ()3.28j) is sim¬ 
ilar, since all the terms in $ are smaller than terms looking like (7(1-1- 

t)~ 2 e , which is bonnded by (7(1 -|- t)“4(l + \x — a^t\ -|- y/t)~ 2 . 

This hnishes the proof of Lemma 13.51 □ 
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